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ABSTRACT: We consider the possibility that small black holes can act as nucleation 
seeds for the decay of a metastable vacuum, focussing particularly on the Higgs po¬ 
tential. Using a thin-wall bubble approximation for the nucleation process, which is 
possible when generic quantum gravity corrections are added to the Higgs potential, we 
show that primordial black holes can stimulate vacuum decay. We demonstrate that 
for suitable parameter ranges, the vacuum decay process dominates over the Hawking 
evaporation process. Finally, we comment on the application of these results to vacuum 
decay seeded by black holes produced in particle collisions. 
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1 Introduction 

We live in a world in which the fundamental properties of matter are manifestly un¬ 
changing on the timescale of our everyday lives. Nevertheless, the recent discovery of 
the Higgs boson [1, 2] raises the possibility that, even within the standard model of 
particle physics, the present vacuum state of the universe may not be stable, but only 
metastable, with another lower energy state at high expectation values of the Higgs 
held [3-7]. In general, this would not conflict with observation because the lifetime 
of the present vacuum would be far longer than the age of the universe. Indeed, the 
possibility that we live in a metastable state was mooted long before the discovery of 
the Higgs [8-16]. 

Investigations of vacuum decay in the context of quantum held theory are usually 
based on the bubble nucleation arguments of Coleman et al. [17-19], (see also [20]) 
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which relate vacuum decay to the random nucleation of critical bubbles of a new vacuum 
or phase. However, in many familiar examples of phase transitions beyond the realm 
of particle physics, the transition is dominated by bubbles which nucleate around fixed 
sites, usually impurities in the medium or imperfections in a containment vessel. It 
is therefore important to investigate whether the metastable Higgs vacuum might be 
ruled out if the seeded nucleation rates for vacuum decay are comparatively large. 

In recent work [21], following earlier work by Hiscock and Berezin [22, 23], we 
looked at the effect of gravitational inhomogeneities acting as seeds of cosmological 
phase transitions in de Sitter space. We found that the decay rates were considerably 
enhanced by the presence of black holes. Following our work, Sasaki and Yeom [24] 
have investigated the unitarity implications of bubble nucleation in Schwarzschild-Anti 
de Sitter spacetimes (see also [25] for a discussion of vacuum stability in the early 
universe). In this paper we extend our previous results [21], to cover all possible 
gravitational nucleation processes, focussing in particular on the nucleation of bubbles 
of Anti de Sitter (AdS) spacetime within a vacuum first reported in [26]. 

We follow the approach of Coleman and de Luccia [19], and assume that the nu¬ 
cleation probability for a bubble of the new phase is given schematically by 

T = Ae~ B i ( 1 . 1 ) 

where B is the action of an imaginary-time solution to the Einstein-Higgs held equa¬ 
tions, or instanton, which approaches the false vacuum at large distances. However, 
unlike Coleman and de Luccia, we consider a spherically symmetric bubble on a black 
hole background. The nucleation process typically requires an instanton that has a 
conical singularity at the black hole horizon. Analogous instantons were considered be¬ 
fore in [27, 28] and fall within the generalised type introduced by Hawking and Turok 
[29, 30]. As in our previous paper, we show that the nucleation probability is well- 
defined. An alternative interpretation of (1.1) and the instanton has been given in 

[31]- 

The vacuum decay process is based on a static black hole, in which a bubble 
nucleates outside the black hole and either completely replaces the black hole with 
a bubble of true vacuum expanding outwards, or nucleates a static bubble leaving a 
remnant black hole surrounded by true vacuum. This latter solution is not stable, and 
small fluctuations will lead it to either expand as with the first situation completing 
the phase transition, or to collapse back inwards leaving the initial state unchanged. 
Of course, this description does not explicitly account for any time dependence of the 
black hole due to Hawking evaporation, however, we can apply the same argument 
as that employed for black hole particle production, namely, we consider only vacuum 
decay precesses which have timescales short compared to the evaporation rate. In other 
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words, we have some confidence in our results when the vacuum decay rate exceeds the 
mass decay rate of the black hole. (The effects of Hawking radiation on tunnelling rates 
have been investigated in [32, 33]). 

We will show that the vacuum decay seeded by black holes greatly exceeds the 
Hawking evaporation rate for particle physics scale bubbles. This clearly has relevance 
for the Higgs potential, which we consider explicitly in §4. A primordial black hole losing 
mass by the Hawking process would decay down to a mass around 10-100 times the 
Planck mass and then seed a vacuum transition. The fact that this has not happened 
therefore means that either the Higgs parameters are not the the relevant range (a 
small region of parameter space for this purely gravitational argument) or there are no 
primordial black holes in the observable universe. 

Since our main application is to the Higgs vacuum, we will first summarize some of 
the features of the Higgs potential relevant to the calculation. As with the phenomeno¬ 
logical explorations of the Higgs potential, we write the potential in terms of an overall 
magnitude of the Higgs, 0, and approximate the potential with an effective coupling 

A e ff, 

y{4>) = ^A eff (</>)</> 4 . (i.2) 

The exact form of A e fr is determined by a renormalisation group computation with the 
parameters and masses measured at low-energy. Two-loop calculations of the running 
coupling [3, 34-36], can be approximated by an expression of the form 

A e fr ~ A* + b ^ln , (1.3) 

where —0.01 < A* < 0, 0.1 M p < </>* < M p and b ~ 10 -4 . The uncertainty on these 
parameter ranges is due mostly to experimental uncertainties in the Higgs mass and the 
top quark mass, however the possibility of negative A e fr approaching the Planck scale is 
quite real. The present-day broken symmetry vacuum may therefore be a metastable 
state, but quantum tunnelling in the Higgs potential determined by the usual Coleman 
de Luccia expressions is very slow, and the lifetime of the false vacuum far exceeds the 
lifetime of the universe. 

The observation of negative A e fr of course assumes no corrections from new physics 
between the TeV scale and the Planck scale. We might expect quantum gravity, or 
other effects will have to be taken into account. On dimensional grounds, we can write 
modifications to the potential of the following form [37-42], 

| A efr(</>)0 4 + ^(<5A)bsm</> 4 + + 8^ 8 ^T + ■ ■ ■ (1-4) 
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Figure 1. The Higgs potential at large values of one of the Higgs held components 0. The 
parameter values for the blue line are A* = —0.001, 0* = 0.5 M p . The black line shows the 
effect of adding a 0 6 term with coefficient Ag = 0.34. 

where (5A)b sm includes corrections from BSM physics, and the polynomial terms rep¬ 
resent unknown physics from the Planck scale. If these coefficients are similar in mag¬ 
nitude, then the small size of A e ff at the Planck scale has the consequence that there 
is an intermediate range of 0 where the potential is determined predominantly by A e fj 
and Ag. 

Quantum tunnelling in a corrected potential has been explored by Branchina et 
ah [39, 40] (see also [43]). They considered potentials with A* ~ —0.1, where the 
potential barrier occurs at 0 <C M p , and they further enhanced the tunnelling rate 
by taking Ag = —2. They claimed a greatly enhanced tunnelling rate, with a lifetime 
much shorter than the age of the universe, however, their discussion did not include 
gravitational interactions. 

The aim of the present paper is to analyse models which appear stable on cos¬ 
mological timescales when using the CDL results, but may become unstable due to 
enhancement of the tunnelling rate by a nucleation seed, which we will take to be a 
microscopic black hole. For this semi-analytic investigation, we consider the nucleation 
with thin-wall bubbles of the true vacuum in an analogous way to Coleman and de 
Luccia. In terms of the Higgs, this thin-wall bubble nucleation requires the potential 
to be relatively shallow at the true vacuum, and this requires a large positive 0 6 term. 
To go beyond this approximation, which allows us to use pure gravitational arguments, 
will require a detailed numerical study that we will present in future work. 

The outline of the paper is as follows. We first review then extend the thin wall 
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instanton method in §2, directly calculating the instanton action in the thin wall limit 
as a function of wall trajectory and black holes masses. In §3 we describe the solutions 
for the instantons and discuss the preferred decay process for a general seed mass black 
hole (including charge). In §4 we apply the results to the case of the Higgs potential, 
and present a full comparative calculation with the decay of the black hole due to 
Hawking radiation. Finally, in §5, we discuss possible extensions to higher dimensions 
and collider black holes. Note we use units in which h — c — 1, and use the reduced 
Planck mass M / = 1/(87 tG). 


2 Thin-wall bubbles 


In this section we describe how to construct a thin wall instanton, along the lines of 
Coleman et al. [17-19], but with the difference that we suppose that an inhomogeneity is 
present. Complementary to our earlier work [21], we apply Israel’s thin wall techniques 
[44] to the bubble wall, and describe the inhomogeneity by a black hole. (In appendix 
B we calculate the instanton action for more general inhomogeneous configurations, 
with the proviso that they be static.) 


2.1 Constructing the instanton 

The physical process of vacuum decay with an inhomogeneity can be represented grav¬ 
itationally by a Euclidean solution with two ‘Schwarzschild’ bulks which have different 
cosmological constants separated by a thin wall with constant tension (for a general 
proof of this result in the context of braneworlds, see [45, 46]). On each side of the wall 
the geometry has the form 


ds 2 = f(r)dr± + 


dr 2 

/Co 


+ r diljj, 


f{r) = 1 - 


2 GM± A ± r 2 


( 2 . 1 ) 


where r± are the different time coordinates on each side of the wall, and the wall, or 
boundary of each bulk, is parametrised by some trajectory r = R( A) (the angular 6 
and 0 coordinates are the same on each side). The Israel junction conditions [44] relate 
the solution inside the bubble with mass M_ and cosmological constant A_, to the 
solution outside the bubble with mass M + and cosmological constant A + . Since the 
bubble exterior is in the false vacuum, we have A + > A_. (A + < A_ was discussed by 
Aguirre and Johnson [47, 48], and the case = 0 has been discussed by Sasaki and 
Yeom [24]). In general, the bubble will follow a time-dependent trajectory representing 
a reflection, or bounce. 
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Following the Israel approach [44], we choose to parametrize the wall trajectory by 
the proper time of a comoving observer, i.e. A is chosen so that 

ffl + y = 1 ( 2 . 2 ) 

and take normal forms that point towards increasing r: 

n± = f± dr± — r dr±, (2.3) 

where dots denote derivatives with respect to A. We also take t± > 0 for orientability 
(see also [24]). In these conventions, the Israel junction conditions are 

f+i+ ~ f-t- = -AnGaR. (2.4) 


The combination of surface tension and Newton’s constant recurs so frequently that 
for clarity we define 

d = 2nGa. (2.5) 


To find solutions to the equations of motion, first note that the junction condition 
(2.4) implies 

/± f ± =(/ ± -^) 1/2 =F^(2, 6 ) 

It is convenient to rewrite this as an equation for R using the explicit forms for f± 


R 2 = 1 


2 A (AA) 2 \ 
d 2 + - + v ' ' 


144cx 2 


) 


R 2 - — [M + 

R 1 


AMAAA (GAM) 2 
24d 2 ) ~ 4i? 4 d 2 ’ 


(2.7) 


where AM = M + — M_ and M = (M + + M_)/2 with similar expressions for A. 

Although this seems to be a more complex system than that considered in [ 21 ], in 
fact it is possible to rescale the variables so that the analysis is very nearly identical to 
that in [21]. To begin with, define 


7 = 


Adi 2 


a 


AA’ ' 1 + Ad 2 i 2 ’ 

and rescale the coordinates to R = aR/'y, t — ar/'y, A 

2 aGM_ (1 — a)aGAM 

h = -+ ---, k 2 

7 a r y z 

gives a Friedman-like equation for R( A): 

AV = 1 _A + Fy_h 


dX 


R 2 


R 


A.y 2 

1+ 3 - 

( 2 . 8 ) 

cxA/ 7 . Then writing 

a 2 GAM 

(2.9) 

2d , y 2 

-U(R) 

( 2 . 10 ) 


together with equations for t± (given in appendix A). These equations with a = 1 are 
precisely the system explored in [21]. The allowed parameter ranges for k\ and k 2 are 
obtained similarly, and discussed in appendix A. 
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2.2 Computing the action 

To compute the action of the bounce, we need to compute the Euclidean action of the 
thin wall instanton: 


Ik — 


167rG 


<Ma 


87 tG 


VhIG 




-2A+) 

87 tG 


167tG 


'M- 


Vg(n_ - 2A_) 


VhK _ + / aVh 


( 2 . 11 ) 


'OM- 




and subtract the action of the background. In this expression, <9A4± refers to the 
boundary induced by the wall - there may also be additional boundary or bulk terms 
required for renormalisation of the action (see below). Note that we have reversed the 
sign of the dM.- normal in the Gibbons-Hawking boundary term so that it agrees with 
the outward pointing normal of the Israel prescription. On each side of the wall in the 
bulk we have 1Z± = 4A±, and the Israel equations give K + — K_ = — 12nGa. 

There are three parts to the computation of the action, _M_, -M + , and W. 

• M~ 

Integrating the bulk term for the ” side of the wall has two contributions, one 
from the cosmological constant in the bulk volume, and a contribution from any conical 
deficit at the black hole event horizon, should one exist. A description of how to deal 
with conical deficits was given in an appendix of [ 21 ], essentially the deficit gives a 
contribution proportional to the horizon area times the deficit angle. Supposing that 
the periodicity of the Euclidean time coordinate, /3, set by the wall solution, may not 
be the same as the natural horizon periodicity, /3_ = 4 7 rr/,,/(! — A_r^), this gives a 
contribution to the action from of: 


Im- — 


/3- — (3 A- 


P- 

A _ 


4 G 


-y / 2 -A(R 3 -rl)dT_ 


1 

AG 


+ A 

AG AG 


A- 

P- 


+ 


'lAzli _ 2 GM_ 


+ 4G / dXR2f - T ~ 


( 2 . 12 ) 


where A- is the area of the black hole horizon in Inserting the value of /3_, and 
taking into account the value of r^, the term in square brackets is identically zero, and 
this contribution to the action does not explicitly depend on the periodicity or indeed 
any conical deficit angle. 

• M+ 

The computation of the action of M.+ is a little more involved, as different regulari- 
sation prescriptions are needed for the different asymptotics of (A)dS or flat spacetime. 
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For Schwarzschild de Sitter, the radial coordinate in the static patch has a finite 
range, and terminates at the cosmological horizon r c , which has a natural periodicity 
P c = - 47 rr c 2 /( 2 GM+ - 2A+r 3 c /3). 


Im+ — 


(Pc ~ P)A C 
4 Gp c 

Ac + JL' 

4 G 4 G 


lp±(r> - R>)dr + 


1 

AG 


A c 

Jc 


3 

2A +rl 


+ 2 GM, 


s? / 


(2.13) 


(2.14) 


where *A C is the area of the cosmological event horizon. Once again, substituting the 
values of f3 c and r c demonstrates that the bracketed term vanishes. For future reference, 
we note the value of the background SDS action at arbitrary periodicity derived in [ 21 ]: 

A c A + 

W - - Jc ~ Jc 

where A+ is the horizon area of the black hole of mass M + . Note that this expression 
is /3— independent as discussed in [21]. 

For Schwarzschild (and Schwarzschild-AdS) the range of the radial coordinate is 
now infinite, and we must perform a renormalization procedure. For Schwarzschild, 
there is no contribution from the bulk integral, and instead we consider an artificial 
boundary at large r o, with a subtracted Gibbons-Hawking term so that flat space has 
zero action [49]. 

i ^=^L n VK{K ~ Ko)= 

Again for future reference, computing the background Euclidean Schwarzschild action 
at arbitrary periodicity (with the same background subtraction prescription) yields 




d\f' + R 2 f + 


(2.15) 


Iesch — ~ + PM + 


(2.16) 


inputting the value of Psch = 87 tGM + . 

For AdS on the other hand, we must subtract off the divergent volume contribution 
[50] by again introducing a fiducial boundary at r 0 , and subtracting a pure AdS integral, 
which must have an adjusted time-periodicity so that the boundary manifolds at r 0 
agree: 


Po = P 


f. 


1/2 


(1 - A + r 0 2 /3)V2 


~ 1 + 


3 GM + 

TqA 


P. 


(2.17) 


Thus 


Im+ — —— 


1 

AG 


dr- 


2A 1 


(r 0 J - R* 


+ AG 


drn 


2A 


+ r 3 
'0 


= ^-4 G 


d\f'R 2 T + 


3 


(2.18) 













i.e. an identical result to the Schwarzschild case (2.15). Computing the background 
Euclidean Schwarzschild-AdS action at arbitrary periodicity we get 

Iesads = — + /3M + (2.19) 

again, the same expression as for Schwarzschild, (2.16). 

• W 

Finally, the contribution to the action from the wall has a particularly simple form 
as the Gibbons-Hawking boundary terms from the wall come in the combination of the 
Israel junction conditions. We therefore obtain 

Iw = ± z~n f 'Sh R + [ °Vh = - f ?-v / h = -^ f dX R (f+f+ - /_r_) (2.20) 

° 7r( - r JdM± Jw JW z J 

having used /+r + — /_r_ = —2 oR. 

Putting all of these results together, we End that the action of the instanton solution 
is 


Ie 


A_ 

4 G + 2 G 


fdM + 

dX [(R - 3 GM + )t+ - (R - 3GM_)r_] + A c 

l _ 4 G 


A+ < 0 
A+ > 0 


( 2 . 21 ) 


Thus the bounce action, given by subtracting the background Schwarzschild/S(A)dS 
action is: 


B = id ~ w + h f iX h 2fl/+ " r2/ 'X f+ ~ (2Rf - - R2 f- ) f -} ( 2 - 22 ) 


This expression is the central result of this section, and is independent of any choice of 
periodicity of Euclidean time, and independent of the choices of cosmological constant 
on each side of the wall. It is in fact even valid when the black hole is charged, as we 
will consider in the next section. 


3 Instanton solutions 

In the previous section we derived the equations of motion for a bubble wall separating 
a region of true vacuum from the false vacuum, and derived the “master expression” 
(2.22) for the instanton action. In this section we discuss general properties of these 
solutions, and demonstrate how the action varies as we change the seed black hole mass 
and the wall tension. Rather than presenting absolute values of the bounce action, it 
proves useful instead to present a comparator to the ‘Coleman de Luccia’ action, by 
which we mean the bounce solution in the absence of any black holes (but with, for 
now, arbitrary cosmological constants). 
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3.1 Coleman de Luccia 

The ‘CDL’ bubble wall satisfies (2.10), (A.l), and (A.2), which are solved by 


R = cos A 

t- = . arctan \/ a 2 — 1 sin A ; 

y/oF^l ( 3 . 1 ) 


(where A G 
Bcdl = ~ 


U = 


a 


\J a 2 — (1 — 2cr 7 ) 2 


yja 2 - (1 - 2 ct 7 ) 2 . \ 

arctan------sm A 

(1 - 2d 7 ) 


|, |] for the full bounce) and the action can be computed analytically as 


2G R<T +- T - > = 


7T 


cry 


G cr(o T l)(ct T 1 — 2cry) 


A + =0' tt£ 2 16 (oT) 4 

G (1-4 a 2 £ 2 ) 2 [ } 


Note that by analytic continuation, these expressions include arbitrary A±, for which 
a < 1 or 1 — 2<77 are possible. In this special case the symmetry of the bubble solution 
has been raised from 0(3) to 0(4), and the result for the tunnelling rate agrees with 
explicitly 0(4) symmetric methods. 


3.2 The general instanton 

The general bubble wall will have a black hole mass term on each side, and a general 
instanton will consist of a bubble trajectory between a minimum and maximum value 
of R. For fixed seed mass, M + , there will be a range of allowed k\ and k 2 (see (A. 5)), 
and a corresponding range of values for the bounce action. By exploring the {k \, k 2 } 
parameter space numerically and plotting the ratio of the bounce action to the CDL 
action, we can build up a qualitative understanding of the preferred instanton for 
vacuum decay. 

For example, if A + = 0, GM + = 7 fc 1 /2a, and GM_ = GM + — 7 fc 2 ( 1 — a) jo?. 
Referring to figure 9, we see there are two possibilities for the range of k 2l which is now 
a horizontal line in the k\ plot: Either the maximal value of k 2 lies on the fcj” branch 
with GM_ = 0, or on the static branch k\{k 2 ). The picture is similar for general A + , 
but the constant GM + lines are now at an angle, and interpolate between the k™ curve 
at negative k 2 and either the k™ line at positive k 2 or the k\{k 2 ) curve. The crossover 
between the two possibilities occurs at M + = Me , given by the algebraic solution to 

?k 

fct(fc 2 ) = -(!-«) (3.3) 

a 

when we have a static bubble with GM_ = 0. In either case, as k 2 drops, GM_ 
increases until the lower limit of k 2 is reached at negative k 2 on the fc” l (/c 2 ) curve. By 
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solving numerically for the wall trajectories we find that the action increases as k 2 
drops. The preferred instanton therefore is the one with the maximally allowed value 
of k 2 consistent with the value of GM + . 

This qualitative picture remains true irrespective of the values of A±: for seed mass 
M + < Mq , the dominant tunnelling process leaves behind a true vacuum region and 
removes the black hole. The tunnelling rate is always faster than the vacuum tunnelling 
rate for these instantons. The bounce action reaches a minimum at M = Me, where the 
bubble is static. For M > Me the dominant tunnelling process is a static bubble with 
a remnant black hole being left behind. As the seed mass increases further, eventually 
the tunnelling rate becomes lower than the vacuum tunnelling rate. Exploring the 
instantons for general A’s, we find that the ratio of B/Bcdl changes very little as the 
A’s vary. In figure 2, we show how this dominant tunneling action varies as the values 
for the cosmological constants are changed. Since the change in B/Bcdl is minimal 
(and Bcdl itself is not varying much), we now restrict our discussion to the A + = 0 
set-up where a = 1 — 2<77, and many of the formulae simplify. 



Figure 2. A plot of the minimum bounce action as M + is varied for at = 0.1, and varying 
values of A_|_ = 6/t 2 , 3 /t 2 , 0, —3/t 2 , A_ = 3/t 2 , 0, —3/t 2 , —6/t 2 as indicated. The ratio of the 
bounce action to the CDL value is plotted, but as A± vary, this value itself changes. For 
at = 0.1, Bcdl = 0.101, 0.117,0.137,0.165 t 2 /L 2 as A + drops from its maximal to minimal 
value considered here. 

Before discussing the general dominant tunneling process, we begin by considering 
the critical instanton where the static bubble tunnels and removes the seed black hole 
altogether. Although (3.3) in general is a complicated algebraic equation, for small at 
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the various parameters can be expanded straightforwardly to give 


he - 



4 

9 


— 3k 2C 


GM C 

i 




From (2.22), the action of a static bounce in general is 


(3.4) 


B = Id ~ 1g = 4?rGM ‘ ~ nG (s) /3 { ^ 3 “ / ‘- 3)2 ' (3 - 5) 

where 

I W GM 

^CP Ml + -±— (3.6) 

although it must be noted that, for the static bubble M_ is a complicated function of 
M + . For the critical bubble, GM_ = 0, hence the critical bounce action is 


B c = AttGMc ~ 


7t£ 2 

~G 




(3.7) 


Thus as at —> 0, the tunnelling action becomes small compared to the CDL action. 

One problem with having a small critical mass is of course that the decay rate due 
to tunnelling may be outstripped by the evaporation rate of the black hole, as we will 
discuss later, however, what this expansion indicates is that the minimal bounce action 
for a particular at can be extremely small, so that even if we are above the critical 
black hole mass, the decay rate can still be significant. 

Having determined that the dominant tunneling process is either the static bubble 
or the GM_ = 0 branch, we can now compute the dominant bounce action either 
by numerically solving the time-dependent bubbles with GM_ = 0, or computing the 
static bubble actions with h = k\. We used a simple mathematica program to calculate 
these exponents, and double checked by a totally numerical computation. The results 
for some sample values of at are presented in figure 3. 

The general bubble solution for GM_ = 0 oscillates between two values Rmax 
and Rmin where the potential U(R) vanishes. This periodic solution in A can only be 
single-valued in if the manifolds on each side have the same time-periodicity as the 
bubble wall solution. In general, this will not be the same as the natural periodicity 
At + = 8nGM + of the Euclidean Schwarzschild solution, hence the need to consider 
general periodicity in the computation of the action in the previous section. For the 
static solution of course, this is not an issue. The values of Rmax/Rmin are well outside 
the black hole horizon radius, and move together as GM + is increased. Eventually, at 
GMc, the two roots of U meet, and the static branch begins. 
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Figure 3. The exponent B for the dominant tunnelling process divided by the appropriate 
vacuum tunnelling value BqdLi for different masses M + of the nucleation seed. The surface 
tension a and AdS radius l enter in the combination at. 

The static branch is the preferred instanton with nonzero GM_, i.e. with a black 
hole remnant, although non-static solutions exist with higher action and remnant mass. 
Initially, the static bubble shrinks with increasing GM + , but remains well outside the 
Schwarzschild radius, however, as we increase GM + further, the bubble becomes con¬ 
strained by the expanding black hole horizon, and becomes stretched just outside Rsch- 
Meanwhile, the remnant black hole mass GM_ increases along the static branch and 
eventually becomes larger than GM + , however, because of the negative cosmological 
constant, the horizon radius, while increasing, does not increase as rapidly as Rsch- 
The static bubble action therefore increases as GM + increases, eventually becoming 
larger than Bcdl (see figure 3). Figure 4 illustrates the behaviour of these mini¬ 
mal/maximal and static values of R as Rsch — 2GM + varies, the remnant horizon 
radius is also shown. 

3.3 Charged black hole instantons 

Finally, before considering the case of the Higgs vacuum in detail, we conclude this sec¬ 
tion by commenting briefly on an obvious generalisation of our instantons to Einstein- 
Yang-Mills-Higgs theory. The combination of Einstein gravity with Yang-Mills and 
Higgs fields admits the possibility of charged black hole solutions [51, 52]. Electrically 
charged black holes can discharge by the emission of charged particles [53], but mag¬ 
netically charged black holes may be the lightest magnetically charged particles in the 
theory, in which case a large mass black hole evaporates towards the extremal limit, 
and the Hawking radiation flux falls to zero. 
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Figure 4. A plot of the variation of the bubble wall radius R as GM + is increased for 
at = 0.25 (chosen to highlight the qualitative features). As at drops, the features of the 
phase diagram remain the same, but ‘squash up’ towards smaller Rsch- The unlabelled red 
line running from corner to corner represents Rsch > the seed black hole horizon radius. 

Magnetically charged black holes may be produced in the early universe [27, 28], 
and form the seeds for vacuum decay of an unstable standard model Higgs held. Un¬ 
charged black holes can easily evaporate before they seed a phase transition, but the 
charged black holes hang around for a longer time making them better candidates for 
vacuum decay nucleation sites. 

An SU( 2) x U( 1) Yang-Mills theory with Higgs held TL in the fundamental SU( 2) 
representation has no hat-space monopole solutions, but it does have Dirac and Yang- 
Mills black-hole monoples. The non-abelian monopoles can be constructed from the 
SU( 2) helds W using an ansatz 

TL = (f)(r)a r TLo, (3.8) 

W = G U 2 — (a a, dd — a e sin 6 d(j >), (3.9) 

r 

where ay, ae and a<p are Pauli matrices projected along the spherical polar co-ordinate 
frame and TLq is constant. (The magnetic charge has been scaled so that an extreme 
black hole has P — M in the absence of a cosmological constant.) 
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For a potential which allows decay from flat space to AdS, there are thin-wall 
bubble solutions with spherical symmetry and constant values of (j) at the appropriate 
minima of the potential. The metric coefficients are 


f- 

f+ 


2GM- | r 2 | G 2 P 2 
/ ^ ~£2 ^ ^2 

2 GM+ G 2 P 2 

- 1 - o- 


(3.10) 

(3.11) 


In this case, the bubble wall carries no magnetic charge. Generalised solutions may 
also be possible in which the interior and exterior have different magnetic charges. 

The action for the bubble solutions is given by the same formula, (2.22), as in 
the uncharged case, though with the appropriate expressions for f±. The plot of the 
dependence of the action on GM + /£ is surprisingly similar to the uncharged case at fixed 
ratio P/M+, with one small modification. The time-dependent tunneling solutions prior 
to the switching on of the static bubbles now do not remove the black hole altogether 
as this would leave a naked singularity. Instead, the bubbles leave behind an extremal 
remnant, M_ = M ext (P ), where 


The static branch meets this time-dependent branch at a critical mass Mcp , where the 
static bubble now has an extremal black hole in its interior. On the static branch, the 
action is, as before, the difference of the areas of the seed and remnant black holes, 
but as the extremal limit is approached, the horizon radius of the remnant black hole 
shrinks only as the root of M + — Mcp , whereas the radius of the seed black hole (which 
is not approaching an extremal limit) depends linearly on M + — Mcp , thus, as we 
increase M + from Mcp , the action actually starts to drop briefly, before the effect of 
the increasing horizon area kicks in causing the usual rising of the bounce action. This 
small dip in the action near the critical point is very hard to see at low P/M + , but for 
larger ratios becomes more visible. The clip is however very minor, and the minimum 
action is well approximated by the value at Mqp- 

From figure 5 we see the dip is most visible at large ratio P/M, however, perhaps 
surprisingly, it is also the case that at large P/M the catalytic effect of the black hole 
is much reduced. We therefore expect that the addition of a monopole charge will 
not particularly assist with vacuum decay, a conclusion largely borne out by the more 
detailed analysis of the next section. 
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Figure 5. The exponent B for the dominant tunnelling process with black-hole monopoles 
of mass M + acting as nucleation seeds. 

4 Application to the Higgs vacuum. 

Up to this point, the vacuum decay process has been described in gravitational terms 
using the surface tension of the wall, c, and the AdS radius of the ‘true’ vacuum, i. 
In this section we will explore vacuum decay in the Higgs model with high energy 
corrections as discussed in the introduction. The key features of the potential relevant 
for quantum tunnelling are the barrier height, the separation between the minima and 
the energy of the true vacuum (TV). These three parameters can be encoded as follows, 

g = (hv/M p , e = -V(<j>rv), C = sup V(4>) (4.1) 

0 <(f)<(f)TV 

Following our previous discussion we shall restrict attention to potentials which 
allow thin-wall bubbles. Although we would expect ( e for a thin wall bubble, 
numerical solutions show that the wall approximation is reasonably accurate even when 
( ~ e, therefore we use this lower bound for The range of Higgs model parameters A*, 
</>* and Ag which allow thin-wall bubbles is set by ( > e > 0, and by the condition that 
the true vacuum lies at large qb. Thin wall bubbles correspond to rather large values 
of Ag, as illustrated in table 1. Roughly speaking, as A* becomes more negative, the 
values of Ag required for thin-wall become larger, similarly as 0* drops, Ag increases. In 
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all cases the pure CDL tunneling action is extremely large (10 6 7 ), but the suppression 
of the critical tunneling action is also large, and increases as A* becomes more negative. 

Table 1. A selection of parameter values for the modified Higgs potential, including the 
AdS radius £, the rescaled surface tension dl and the critical mass Me for optimal nucleation 
seeded by a black hole. These parameters lie along the bottom edge of the parameter ranges 
for thin-wall bubbles. The vacuum tunnelling exponent Bcdl is around 4e+06 in each of 
these examples. 


Aj): 

4>*/M p 

A6 

9 

£/L p 

a£ 

Mc/M p 

Be/B C dl 

-0.005 

1 

500 

0.00146 

3.17e+8 

0.00045 

3.5 

1.2e-6 

-0.005 

0.5 

2e+03 

0.00073 

1.27e+9 

0.00023 

1.8 

2.9e-7 

-0.007 

2 

1.98e+03 

0.0008 

9.33e+8 

0.00024 

1.5 

3.2e-7 

-0.007 

1 

7.93e+03 

0.0004 

3.79e+9 

0.00012 

0.8 

8.2e-8 

-0.007 

0.75 

1.41e+04 

0.0003 

6.76e+9 

9.1e-05 

0.61 

4.7e-8 

-0.007 

0.5 

3.17e+04 

0.0002 

1.51e+10 

6e-05 

0.39 

2e-8 

-0.008 

1 

27e+03 

0.00022 

1.18e+10 

6.9e-05 

0.46 

2.7e-8 

-0.008 

3 

3e+03 

0.00067 

1.31e+9 

0.00021 

1.4 

2.4e-7 

-0.009 

1 

85e+03 

0.00013 

3.43e+10 

4.1e-05 

0.28 

9.4e-9 

-0.01 

2 

63e+03 

0.00016 

2.55e+10 

5.4e-05 

0.49 

1.7e-8 


Following Coleman and De Luccia, we can express the surface tension of the bubble 
wall in terms of the potential. In order to extend the result to moderate values of £/e, 
we compute the tension using the integral 

d(p(2V) 1/2 ~ KgM p ( 1/2 , (4.2) 

where the upper limit of the integral is at V((j> i) = 0. The constant k depends on the 
details of the potential, but since (j)\ < gM p and V < £, it is subject to the constraint 
k < a/2. The AdS radius i is related to the vacuum energy density by 

l = VSM p e~ 1/2 . (4.3) 



The back-reaction parameter a£ is therefore 


d£ = -r^a£= ^ 


4 Ml 


1/2 


(4.4) 


Note that di < 1/2 puts an upper bound on g. The CDL tunnelling exponent Bcdl 
given in (3.2) is 

r 2^4 nn.. 4_2 / a \2 


Bcdl — 


Th ru 4 


2e 3 (l - 4 d 2 £ 2 f 


27 (4 . 5) 
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The large size of Bqdl in the parameter range covered by table 1 guarantees a tunnelling 
lifetime longer than the age of the universe (for unseeded nucleation). 

When the vacuum decay is seeded by a black hole, the most rapid decay process 
occurs for a seed mass M + = Me given in (3.4), 


M c 


128 i 
~27G 


16 


,crn = 


-'KK 




(4.6) 


The corresponding exponent in the nucleation rate is Be = 0.5 (Me / M p ) 2 . Some values 
for Me are shown in table 1. If Me M p , then the exponent Be is large and the 
seeded decay rate becomes vanishingly small. On the other hand, if Me < M p then 
even if we have a seed mass M + Me, we can still get significant suppression of 
the bounce action while remaining well above the Planck mass from tunneling on the 
static branch, as we see from figure 3. Our strategy therefore is to explore the decay 
seeded by a small mass black hole via the static instanton, which (for convenience) we 
determine numerically. 

A brief consideration of the dependence of the bounce action on M + shows that 
we are exploring seeded tunnelling for very light or primordial black holes [54], with 
temperatures well above that of the CMB. We must therefore check that the black hole 
can seed the false vacuum decay before it disappears through Hawking radiation. The 
vacuum decay rate T^, (1.1), contains not only the exponential of the bounce action, 
but also a pre-factor, A. According to Callan and Coleman [18], this pre-factor is made 
up of a factor of (H/27T) 1 / 2 for each translational zero mode of the instanton and a 
determinant factor. In our case, there will be a single zero mode representing the time 
translation symmetry, and rather than evaluate the determinant factor, we use the 
inverse horizon timescale as a rough estimate (GM + )^ 1 , giving 



(4.7) 


The black hole emits Hawking radiation at a rate depending on fundamental particle 
masses and spins. The total decay rate for a subset of the standard model was evaluated 
by Page, [55]. If we set T H = M/M, then 


r„ ss 3.6 x 10- 4 (G 2 AfJ)-‘ 


(4.8) 


The branching ratio of the tunnelling rate to the evaporation rate for uncharged black 
holes is therefore 


Td 

r H 


44 


Ml 


M 2 


B 1 / 2 e~ B . 


(4.9) 
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Figure 6. The branching ratio of the false vacuum nucleation rate to the Hawking evapora¬ 
tion rate as a function of the seed mass for a selection of Higgs models from table 1. The first 
plot shows the branching ratio for 0* = M p with the labelled values of A*, and the second for 
A* = —0.007 for the labelled values of 0*. The black hole starts out with a mass beyond the 
right-hand side of the plot and the mass decreases by Hawking evaporation. At some point, 
the vacuum decay rate becomes larger than the Hawking evaporation rate. 

From this expression we can see roughly how the branching ratio will depend on 
M .|_, even though B is, in principle, a complex function of M + . The static bubble is 
the difference in areas of the seed and remnant black holes, which we can deduce from 
figure 4 to be roughly linear (there is actually a slightly stronger dependence on M + , 
however, what is important is that it is not quadratic), whereas the prefactor is (again, 
roughly) M+ 2 ; we therefore expect the plot to be strongly exponentially suppressed at 
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large M + , but rising as M + falls to a maximum around M + /M p = 0( 1), then dropping 
again below M p . The actual value of the maximum will depend on the details of how 
B depends on M + , which requires a full calculation. 

The branching ratio is plotted as a function of the seed mass M + for two indicative 
sets of parameters taken from table 1 in figure 6 in order to illustrate the dependence on 
the parameters in the Higgs potential (or correspondingly on a and €}. The branching 
ratio is shown at fixed 0* with varying A* and vice versa. The overall picture is that 
for lower a and higher i (or more negative A* / higher 0*) the branching ratio is larger, 
and is consistently higher than unity over a larger range. While Hawking evaporation 
always dominates at large M + , the effect of Hawking radiation is that the black hole 
loses mass, hence driving it towards increasing branching ratio. A black hole produced 
in the early universe, for example, starts out with a mass well beyond the right-hand 
side of the plots, but at some point after evaporation, the vacuum decay rate becomes 
larger than the Hawking evaporation rate and the black hole seeds the transition to a 
new vacuum. This can occur for seed masses well above the Planck mass, where we 
have some confidence in the validity of the vacuum decay calculation. The timescales 
for Hawking evaporation and vacuum decay will both be less than roughly a million 
Planck times. 

Finally, for the case of a monopole charged black hole, we might expect the branch¬ 
ing ratio to be larger due to their reduced Hawking radiation rate: The Hawking flux 
is proportional to A+T+, where A+ is the event horizon, and T + is the Hawking tem¬ 
perature 


T + 


1 4A 
8itGM + (1 +A) 2 ’ 


(4.10) 


setting A 2 = 1 — P 2 jM\. The evaporation rate is now 


r 


H 


3.6 x 10 


r3\ —1 


64A 4 

(1 + A ) 6 


(4.11) 


The branching ratio T^/Th can now be re-computed using this evaporation rate and 
the vacuum decay rates from the previous section. The result is shown in figure 7. 


5 Conclusions 

The main aim of this paper was to demonstrate that black holes can massively speed 
up the rate of decay of a metastable Higgs vacuum by acting as nucleation sites. We 
have shown that is the case, subject to the limitations of the analysis, some of which 
we shall now address. 
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Figure 7. The branching ratio of the false vacuum nucleation rate to the Hawking evap¬ 
oration rate for a monopole charged black hole with P = 5 M p , shown as a function of the 
seed mass for a selection of Higgs models from table 1. The plot for the uncharged black hole 
(P = 0) is repeated for comparison. As before, A* = —0.007 for the labelled values of </>*. 

We have used a Higgs model which was based loosely on the two-loop effective po¬ 
tential, with an extra term motivated by high-energy physics, possibly quantum gravity. 
This has enabled us to simplify the analysis by employing a thin-wall approximation, 
which is valid in part of the parameter space. The thin-wall approximation shows that 
when the seed mass is above some (calculable) critical mass, then the seeded nucleation 
proceeds via a static bubble solution. This gives a good starting point for an analysis 
of the thick wall bubble nucleation, which is far simpler for static than for non-static 
field configurations. We will show in a companion paper that the thick-wall bubble 
solutions extend the allowed range of parameters to small values of our new coefficient 
Ag, which was rather large here, to the limit where Ag = 0, and the potential becomes 
simply the two-loop Higgs potential from the standard model. 

The seeded nucleation calculation presented here requires a black hole, and for 
large tunneling enhancement, this is expected to be a primordial black hole which is 
evaporating and nearing the end of its life. There is of course another situation in 
which small black hole might occur. A possible alternative solution to the hierarchy 
problem has been to consider large extra dimensions, [56-59]. In these models, our 
universe is presumed to be a “brane” living within higher dimensions on which standard 
model physics is confined. The higher dimensional Planck scale is not hierarchically 
large, but instead our 4D Planck scale, which is derived via an integration over these 
extra dimensions, gains its leverage via the “large” internal volume. In such scenarios 
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black holes can be created in particle collisions [60, 61], leading to considerable interest 
in the possibility of black holes being produced at the LHC (for a recent review see 
[62]). There are no known exact solutions for these black hole plus brane systems, and 
instead the black hole is usually considered to be approximately a higher dimensional 
Schwarzschild or Myers-Perry [63] solution (see [64, 65] for reviews on the issues and 
properties of brane black holes). For the Randall-Sundrum braneworld models, where 
the extra dimension is strongly warped, one could also consider “brane only” solutions, 
such as the tidal black hole [66]. 

Calculating the vacuum decay rates for these systems would be challenging, to say 
the least, not only because of the lack of a true higher dimensional black hole solution, 
but also because an instanton presumably would have to have a different vacuum only 
on the brane, and not in the bulk (although the braneworld equivalent of the CDL 
instantons were constructed in [46, 68]). However, some features of our calculation 
should be present. The tidal black holes, for example, resemble black-hole monopoles, 
but with negative square monopole charge P 2 * . The bubble solutions for tidal holes are 
simple generalisations of the ones we have already discussed. 

More concretely, in analogy with the use of the Myers-Perry geometries to approx¬ 
imate the small collider black holes, we could consider a higher dimensional analogue 
of our instanton. This would be a flawed model as it does not fit easily with a Higgs 
held which is confined to the brane, but it provides a concrete calculational method 
presumably as reasonable as the use of higher dimensional black hole geometries to 
approximate the brane black hole. Appendix C gives the details of the computation 
with the higher dimensional instanton, which leads to a dimension dependent branching 
ratio 


r H 


550 


M + r + 
(D - 3) 4 



550 y/B^e- BD 


IQtiGdM, 


D—2 


(D-3) 4 l(D- 2)A d _ 2 J 


l/P-3) 


(5.1) 


where Bn is given (implicitly) in (C.6). Defining the higher dimensional reduced Planck 
mass as 1 

<- 2 = s^- d (5 - 2 > 

we can track the branching ratio as a function of Ad + /Mn and its dependence on D. 
To illustrate this dimensional dependence, we chose test-case values of ai = 0.01 and 
£ = 0.1, and plotted Td/Th in figure 8. 

1 Note that in the literature, see e.g. [64], the non-reduced Planck mass is often used. Due to the 

dimension dependence of the Planck mass this will introduce various dimension dependent renormal¬ 

isation factors between our expressions and those assumed there. Although these are of order unity, 

they do have some impact. 
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Figure 8. The dependence of the branching ratio on the dimensionality of spacetime, D. 
Here, the D —dimensional Planck mass is fixed, i.e. 8nGr) always has the same value. 

The figure shows that as the number of extra dimensions increases, the branching 
ratio decreases, however the exact normalisation of the plot will depend on the confi¬ 
dence of translating modifications of the Higgs potential to the new Planck scale and 
variables. Given the crudeness of this particular model, we leave this, and possible 
refinements of the decay modelling to future work. 

In conclusion, we have shown that the lifetime of our universe in a metastable Higgs 
phase is crucially dependent on the absence of any nucleation seeds, and a primordial 
black hole could drastically reduce the time it takes to decay onto a different ‘standard 
model’. Instability of the standard model is therefore more problematic than was 
hitherto supposed. Further exploration of the parameter space, using a wider class of 
bubble nucleation scenarios, should give us the range of Higgs parameters which lead 
to a long-lived standard model in the presence of black holes. 
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A The limits on Aq, /t2 

i 2 

The values of k\ and k 2 are limited by requiring existence of a solution to R +U(R ) = 0, 
where U is defined in (2.10), requiring positivity of the black hole masses, and positivity 
of the arrows of time on each side of the wall: 


where 


/dw = -| + -d- 2 ^) 

d\ R 2 a 
df_ _ k 2 R 
~ d\~ R 2 a 


„ k\ 2 ko . , R 2 r 9 , 91 

/• - 1 - ---(«-(!- 2CT7))-2 [a - (1 - 2cry)-] 

R aR O' 

„ k\ 2ko, . R 2 . 9n 

/- = 1 —-- H-g(l - a) H-y(l - a ) . 

R aR a" 1 


(A.l) 

(A.2) 


(A.3) 


The first requirement is algebraically identical to the constraint discussed in [21] (al¬ 
though the expression given there in the appendix was not correct). Simultaneously 
requiring U = U' = 0 and eliminating R* gives an upper limit k\ < k*, the correct 
expression for which is: 


k* = - 
1 9 


1 + 81Jfe; 


-1 - 5(27k 2 ) 2 + 


1 + 5(27 k 2 f 


(27k 2 

K2) t- 

(27 k 2 ) A , 27 k 2 


+ 


27h 


(4+(27fc 2 ) 2 ) 3/2 ) 


1/3 


+ 


(4 + (27k 2 f) 


2 \ 3/2 


1/3' 


1/2 




2ko 


(A.4) 


To get lower limits for Aq we have to consider positivity of the black hole masses 
and the arrows of time on each side of the wall. These now depend on the sign and 
magnitude of the cosmological constants and are different to [21]. First, note the 
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relation between the physical quantities and the parameters: 

GM_ = (k x - 2 ^ ~ a h 2 

2a V « 

GM + = ^ + (a - 1 + 2cry) 
a- 


A_ = 3 


2 a 

a 


2 -l 


(A.5) 


7 


A + = A_ — 12a 2 + 12- = 4 (« 2 - (1 - 2h7) 2 ) 

ry ry* ' ' 


thus 


GM_ > 0 

GMx > 0 


fci > 2(1 — a) — 
a 

fci > 2(1 — a + 2 (t 7 ) — 
a 


(A.6) 


Secondly, the requirement of positivity of f±, of which the constraint on r + is the 
stronger: 

R 

(A.7) 


ko R, 

2 + _ l_2a 7 >0 
i ? 2 a 


saturated by = ak 2 /{2a 7 — 1). Clearly, if 2(77 > 1, we must have k 2 > 0, and [/ 
must be positive with positive gradient at R + . A brief manipulation of U' > 0 yields 

« 2 ki + 2ko a „ . . 

(2d7 - l) 2 _ 2fc 2 (2d7 - 1) ~ ° ^ 

From this, we see that A + > 0, 2a 7 < 1 + a/2, and k\ is bounded below by /c / 1 = 
2(1 — a)k 2 /a from (A. 6 ). 

Now consider 2(77 <1/2, for which k 2 < 0 is allowed. A similar argument to [21] 
gives 


ki>k™ = 


ak 2 


2(77 — 1 


1/3 


+ 


k 2 (a - 1 + 2 a 7) 5 


(A.9) 


(1 — 2(77) a 

for k 2 < 0, and once again, k™ = 2(1 — a)k 2 /a for k 2 > 0. Note that neither of these 
bounds requires a particular sign for A_. 

Where the sign of A_ does make a difference is in the range of allowed k 2 . The 
upper limit on k 2 is determined by M_ = 0 for the static bounce, i.e. 


akl[k 2 \ = 2(1 — a)k 2 


(A. 10) 


and a lower bound on k 2 is determined when the range of k\ for k 2 < 0 closes off at 
negative k 2 , which occurs when U(R + ) = U'(R + ) = 0: 


R% - 3 R% + 3 k\ = R\ 


^k\ 


(1 — 2ay) 2 


A, = 0 


(A.ll) 
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Figure 9. The allowed ranges for the parameters k\ and Ay with A+ = 0. The upper 
bound on Ay corresponds to stationary wall solutions. The lower bound on Ay when Ay > 0 
corresponds to vanishing remnant mass M_ = 0. The two limits intersect at a point (Aye, Aye), 
which depends on the ai and approaches (0,4/27) as al —> 0. The strip of allowed k\ continues 
to the left as 17c21 1 ' /3 - 


which is clearly inconsistent for A + < 0. For A + > 0, the range closes off at 

, _ a 2 (2d-7-l) _a 2 (2cr7-l) , AloA 

2 ’ mm 3 v / 3[« 2 - (1 - 2 d 7 ) 2 ] 3 / 2 7 3 Af 1 j 

Note that for A + < 0, the range of Ay, while initially narrowing as Ay becomes nega¬ 
tive, eventually opens out, as the linear term in (A.9) becomes dominant. Thus large 
AdS black holes can tunnel to an even larger AdS black hole with a more negative 
cosmological constant. 

In this paper, we mostly focus on A + = 0, for which a = 1 — 2a r y, and many of the 
expressions above simplify: 


K 


(-k 2 ) 1/3 Ay < 0 

^Ay k 2 > o 


(A.13) 


and the range of k 2 is plotted in figure 9. Note that unlike the figure in [21], where the 
lower limit for Ay was dependent on a for negative but not positive Ay, in this case it is 
the lower limits for positive Ay and not negative that depend on a. 


B Alternative bounce action calculation 

In this appendix we present an alternative derivation of the general expression for the 
bounce action using the Hamiltonian approach presented in [21], and extend the re- 


26 












suit in some cases beyond the thin-wall limit. We will evaluate the Euclidean action 
for gravity plus a scalar held, with Lagrangian C rn , in an asymptotically AdS or hat 
spacetime. As in §2.2, we take the action on a sequence A4 r of manifolds with bound¬ 
ary dAi r at ro, and subtract the action of a similar sequence M. r of Euclidean AdS 
manifolds with the same boundary dM. r . 

It is instructive to hrst consider the case where the Euclidean spacetime A4 is 
perfectly regular, with no conical singularities, and has a Killing vector, d T . We perform 
a foliation of A4 r with a family of non-intersecting surfaces E r (assuming the global 
topology permits), with 0 < t < /3. The canonical decomposition of such foliations has 
been investigated by Hawking and Horowitz, [67]. In order to decompose the action, 
we use their identity 

77 = 3 77 - K 2 + K 2 ab - 2 V a (u a V b u b ) + 2V b (u a V a u b ), (B.l) 


where the vector is normal to E T and 3 77 is the Ricci-curvature of E r . The action 
therefore splits into bulk and boundary parts, 

W = - f dr [ ( 3 ^ - K 2 + K ab + 16tt GC m ) ^ 

J 0 St 

+ [ n b u a W a u b Vh. (B.2) 

JdMr 


The bulk term expressed in terms of canonical momenta and n becomes 


1 

167tG 


dr 


(d T 'y i j'K lj + d T cf) 7r — N% — N l Tii) 


(B.3) 


where 77 and 77, are the Hamiltonian and momentum constraints respectively. The 
held equations imply that 77 = 77j = 0, and furthermore the symmetry implies d T (j) = 
dr'jij = 0. Therefore only the boundary term in (B.2) survives. 

To evaluate the boundary term, we use the fact that the metric is static and 
asymptotically AdS, therefore at large r we have 


ds 2 = f(r)dr 2 + f(r ) 1 dr 2 + r 2 dQ 2 


(B.4) 


where / = 1 — 2GM/r — Ar 2 /3. For this metric n b u a V a u b = f 1 / 2 / / /2, and subtracting 
the Euclidean AdS action from (B.2) we arrive at the following expression 
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where f3 0 and f' Q are the time-period and metric function of AdS space, and using (2.17) 
the Euclidean action (B.5) becomes, 

I = PM. (B. 6 ) 


If there is a horizon at r = r/,, then by properly treating the conical singularity as in 
[ 21 ], we get an additional area term contribution to the action 

i = m - (b.7) 

This result generalises a previous result of Hawking and Horowitz, who found the same 
formula for the Euclidean action of static Einstein-matter solutions with A = 0 and 
no conical singularities [67]. The A —>■ 0 case can also be obtained using the Gibbons- 
Hawking subtraction procedure described in §2.2. The expression (B.7) is of course the 
same as (2.16, 2.19), but we have not assumed anywhere for the static case that the 
bubble is in the thin-wall limit. 

Solutions with a moving bubble wall (r(A),A(A)) break the time-translation sym¬ 
metry of the full space-time, but the canonical method can still be used if the wall is 
thin and the geometries on both sides of the bubble wall, Al±, individually possess the 
Killing vector, d T . Along with the contributions from and the wall, W, we will 
have an additional contribution from the conical singularity which can be dealt with 
by the methods of [21], The action splits into contributions from each of these parts 

I Mr = /- + /+ + Iw ~ (B. 8 ) 

where AL is the area of the black hole horizon in the interior of the bubble. In the thin 
wall limit lyy = j w <J, and 

4 =-ii IgL ± k ^-L Cm(s ’ ^ + 8 kL x± A (B - 9) 


Hence, performing the same decomposition as in the static case we can cancel the bulk 
contributions and are left only with the boundary terms 


I± = 


8 v iG 


K+Vh + 


'w 


87 tG 


n± b u a V a u b Vh, 


(B.10) 


'8M± 


where dM. + now includes the wall and the large distance boundary dM. r . Using Israel’s 
junction condition to relate the extrinsic curvatures on each side of the wall to the 
tension, and inserting the normal to the wall (2.3), n b u a S7 a u b = r/ _ 1 // 2 / / /2, we reach 
our hnal result 

/ = f)M + - ~ l j w ^‘- jPg L - f'M ^ ( B ' n > 
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The expression for the bounce action B which governs the decay rate is obtained from 
/ by subtracting the background action without the bubble, R: 


J 0 = (3M + - 



(B.12) 


Therefore the tunnelling rate is determined by 


p — _ I 

~ 4G 4 G 2 


cr 


Vh 




167tG 


'W 


(f'+T + -/lf_) Vh, 


(B. 13) 


which is identical to (2.22), after using the relation f+f + — /_f_ = —2 crR for the wall 
integral. 


C Higher dimensional instantons 


Here we briefly outline the simple higher dimensional instanton model. Note this is only 
a crude calculation meant to give an indication of the effect of large extra dimensions, 
and is problematic as an approximation to an actual Higgs tunneling event due to the 
lack of an exact braneworld black hole solution, as well as the fact that the Higgs should 
be confined to the brane, hence the instanton should instead involve only a shift in the 
cosmological constant on the brane. We take a higher dimensional black hole solution 
with different masses and cosmological constant on each side of the wall, and compute 
how the branching ratio changes with dimension. 

The equations of motion have the same schematic form as (2.6) 

R = a R -/+ W - 2R2 (C.l) 

but with a = 4nGa/(D — 2) [46, 68 ], and / now the higher dimensional Schwarzschild 
potential: 

2A r 2 16ttG d M 

~ (D — 1 ){D — 2) — (D — 2 )A D _ 2 r D ~ 3 

where Gd is the higher dimensional Newton constant and Ap -2 = 2 it /T[^A]. 
Defining £, 7 and a in a similar fashion: 


(C.2) 


and 


(D-l)(D-2) 

2AA 

k\ 

h 


7 = 


4 a£ 2 


1 + 4a 2 £ 2 


cr = 1 + 


2A 7 2 


(D-l)(D-2) 


(C.3) 


167tG 
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(D — 2)A d _ 2 \7 
I67 tGd 


d -3 


. AM 


(D - 2 )A d _ 2 V 7 


a \ D ~ 2 AM 


(C.4) 


4cr 
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Then setting R = aR/'f, A = a A /7 gives the equation of motion 



X _ + -^2 _ k '2 

RD -3 i?2(D-2) 


(C.5) 


which is of the same form as (2.10) albeit with different exponents of R. We can use the 
same procedure as before to find the static solution and the dynamical bubble which 
removes the black hole. The static action (which is what is needed for the branching 
ratio) is, as before, the difference in horizon areas: 

Bo = (r ?- 2 - r D - 2 ) (C. 6 ) 

where r± are determined numerically, and the corresponding tunneling rate is 


r D = \l — 


B n e~ BD 


(C.7) 


27r r + 

Meanwhile the Hawking temperature of the higher dimensional black hole is [63, 64] 

(D-3) 


Th = 


47rr 4 


(C. 8 ) 


To estimate the decay rate of the black hole, we will assume that the main channel is 
due to emission of particles on the brane [69, 70], leading to 


H 


3.6 x 10 


_ 4 64t tG d (D- 3 ) 4 
Ad-2T?~ 1 (D - 2) 


(C.9) 


hence a branching ratio of 


« 550 M + r + 

T h (D- 3)4 


\/ B*( 


-B , 


(C.10) 
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